Mon. Not. R. Astron. Soc. 000, [T]-?? (2006) Printed 5 February 2008 (MN IATeX style file v2.2) 



Excitation and visibility of slow modes in rotating B-type 
stars 



o 
o 

> 

o 



> 



o 

o 
U 



W. A. Dziembowski^'^*, J. Daszynska-Daszkiewicz^ and A. A. Pamyatnykh^'^ 

^ Warsaw University Observatory, Al. Ujazdowskie 4, 00-478 Warsaw, Poland 
^Copernicus Astronomical Center, Bartycka 18, 00-716 Warsaw, Poland 

^ Instytut Astronomiczny, Uniwersytet Wroclawski, ul. Kopernika 11, 51-622 Wroclaw, Poland 
'^Institute of Astronomy, Russian Academy of Sciences, Pyatnitskaya Str. 4-8, 109017 Moscow, Russia 



Accepted 2006 October 2; Received 2006 September 29; in original form 2006 July 7 



X 



ABSTRACT 

We use the traditional approximation to describe oscillations with frequencies com- 
parable to the angular rotation rate. Validity of this approximation in application to 
main-sequence B stars is discussed. Numerical results regarding mode stability and 
visibility are presented for a model of the Be star HD 163868. For this object, Walker 
et al. (2005) detected a record number of mode frequencies using data from the small 
space telescope MOST. Our interpretation of these data differs from that of Walker 
et al. In particular, we interpret peaks in the lowest frequency range as retrograde g 
modes. We find instability in a large number of modes that remain undetectable be- 
cause of unfavourable aspect and/or effect of cancellation. There is no clear preference 
to excitation of prograde modes. 
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1 INTRODUCTION 

Results of linear stability calculations suggest that excita- 
tion of high-order g modes in main-sequence B-type stars 
must be a common phenomenon. The instability is typically 
robust. There is a significant excess of driving over damping 
and many modes are simultaneously unstable in individual 
stellar models. With the standard Population I composition, 
the instability is found in models of nearly all B-type stars. 
Yet, the observed pulsation associated with the slow modes 
is unspectacular. Individual modes never attain high ampli- 
tudes. Only in a small fraction of B stars do we have a firm 
evidence for the presence of oscillations and typically only 
few modes are seen. One possible reason why, despite strong 
instability, we never see high amplitude pulsation of this 
type is a collective saturation by a large number of modes 
which are difficult to detect. Recent report by Walker et 
al. (2005) on the detection of about sixty periods in the data 
on HD 163868 obtained from MOST may support this idea. 

The fact that HD 163868 is a rapidly rotating Be type 
object makes the detection of many oscillation modes partic- 
ularly interesting. The cause of activity observed in stars of 
this type is still not understood. Seismic constraints derived 
from mode frequencies on the internal structure would be of 
great value in this context. Furthermore, it is possible that 
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oscillations play an active role in these objects by inducing 
angular momentum transport, as first suggested by Osaki 
(1986) and Penrod (1986) and more recently, in the context 
of opacity-driven modes, by Lee (2006). Let us add that an 
asymmetry in excitation of prograde and retrograde modes 
may induce a net helicity of motion and thus an action of 
magnetic dynamo. 

To extract seismic information, as well as to assess pos- 
sible active role of oscillation, we have to identify modes 
excited in the star. The task of individual mode identifica- 
tion in an object like HD 163868 is formidable and will be 
attempted in this work. Like Walker et al. (2005), we will 
only try to determine angular degrees and azimuthal orders 
of the three groups of modes which are centered on periods 
of 8 days, and 14 and 7 hours. These authors, however, have 
not found unstable modes that could explain the longest 
period group. 

We will argue that the observed modes represent only a 
small subset of the excited modes. Other modes having sim- 
ilar rms amplitudes may be not detectable because of can- 
cellation and/or unfavourable aspect. Our analysis of sta- 
bility and visibility is based on linear non-adiabatic calcu- 
lations with the effects of rotation treated in the traditional 
approximation. The same approximations were adopted by 
Townsend (2005a, b) in his extensive survey of stability of 
slow modes in models of rotating main-sequence B-stars. 
The slow modes include high-order g modes with frequen- 
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cies comparable to the angular rate of rotation and certain r 
modes (called mixed gravity- Rossby modes), which at suffi- 
ciently high rotation rate become propagatory in the radia- 
tive envelopes. Both types are candidates for identification 
of modes excited in HD 163868. 

Our identification differs from that of Walker et 
al. (2005). One may suspect that the difference should be 
blamed to our use of the traditional approximation, which 
they avoid. However, as we explain in the next section, the 
approximation should be valid in the application to, at least, 
certain modes in HD 163868. 

In Section 3, we present results on stability of slow 
modes in a model of HD 163868 which is very similar to the 
model used by Walker ct al. (2005) . Wc review properties of 
unstable modes in this and similar models and compare our 
results with those of our predecessors. In Section 4 we quan- 
tify the dependence of observable amplitude on the mode 
geometry and the aspect angle for the most unstable low- 
degree modes. Comparison with observations is presented in 
Section 5. 



2 ON THE TRADITIONAL APPROXIMATION 

There is a number of recent papers where the traditional ap- 
proximation is used to describe properties of low frequency 
modes in rotating stars (e.g. Bildstein, Ushomirsky & Cut- 
ler 1996; Lee & Saio 1997; Townscnd 2003a; Savonijc 2005; 
Townsend 2005a,b). Readers arc referred, in particular, to 
Townsend (2003a) for details of derivation of simplified equa- 
tions for adiabatic oscillations known as Laplace's Tidal 
Equations, which employ the traditional approximation. The 
last three papers extend the approximation to the case of 
non-adiabatic oscillation. Here, we focus on the validity of 
approximation and its extension to the case of non-uniform 
rotation and we will quote only the most essential equations. 

The generally adopted, initial simplifications of the lin- 
earized Euler equations. 



^■Vp-I- V*' = 0, 



(1) 



where all notation is standard, are (i) slow and (ii) uniform 
rotation, (iii) oscillation frequency, of the same order as 
the angular velocity of rotation, f2, (iv) neglect of 3>' (the 
Cowling approximation), (v) adiabaticity of perturbations. 
Here we adopt (i) and (iii), that is, we assume 



' s e < 1, 

9 9 



(2) 



but not (ii) as we allow the r dependence in Q. Indeed, the 
traditional approxiniatiou applies to this more general case 
in the essentially unchanged form. We begin assuming (v) 
but later we will consider non-adiabatic effects, which are in 
fact of main interest for us in this paper. 

Inequality (2) implies that the effects of centrifugal dis- 
tortion may be neglected which together with (iv) reduces 
Eq. (1) to 



dvV 1„ , g' 

dt ) Q Q 



0. 



(3) 



Since we allow non-uniform rotation, the constant eigenfre- 
quency, Wobs, must be calculated in the inertial system, while 



u) - the frequency calculated in the local corotating system 
- may vary with r. The acceleration in the inertial system 
expressed in terms of the displacement vector. 



i = rC(r, 7?) exp[i(mvJ - ojobst)], 



(4) 



is given by 
{ dv\' 



dt J 



dn^ 



= —oj ^—2iD,ojez'x^+{erSm'&+e^cos'&)— — rsini?^,., 

(5) 

where uj{r) = u!ohs—'mil{r). Note that with the adopted form 
of the time dependence, the modes of azimuthal order m > 
are prograde and those with m < arc retrograde. This 
is opposite to what has been adopted in the cited papers 
employing traditional approximation but it is in agreement 
with the (7e§Tne Resolution (see Appendix to Christensen- 
Dalsagaard & Dziembowski 2000). 
Using the adiabaticity condition 



Q Vi p r 



(6) 



where we denoted 
A = - 



dlng _ 
dlnr 



and Vg 



ggr_ 



and introducing one additional eigenfunction, Cp{r,'&), 
through the following identity 

p' = grgCpir, ■&) exp[i(mip - Wobat)], (7) 

we get from Eq. (3) 

^N^ - -h sin^ Cr + Sifiw sin i^Cv + 

f ^|: + t^-^-i)Cp = 0, (8) 



dr 
and 



r sin cos t?Cr -u!^Ci> + 2iOw cos t?Cv + - ^ = 0, (9) 

r ov 



2fla; sin -I- 2Q.u) cos -OC,^ — \oj C,^ 



mg 



■Cp = 0. (10) 



r sin t? 

In Eq. (8) we use standard notation A*' for the Brunt- Vaisala 
frequency and U for the logarithmic derivative of the frac- 
tional mass. The continuity equation combined with Eq. (6) 
yields 

\ im 



or 



(11) 

which closes the system of equations for the eigenvalues Wobs 
and the associated eigenf unctions i^r, C^, Cv; ''■nd Cp- 

In the traditional approximation all terms arising from 
the Coriolis force, except those in Eqs. (9) and (10) which 
contain horizontal components of are dropped out. A sim- 
ple justification, based on the local wave approximation, was 
given by Lee and Saio (1997). The approximation is valid if 



N = 



(12) 



Earlier, the same condition was obtained by Dziembowski 
& Kosovichev (1987), who did not rely on the wave ap- 
proximation and did not assume uniform rotation. Condi- 
tions (2) and (12) justify the neglect of those Coriolis terms 
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and the term containing the derivative of fl, both in the 
outer evanescent zone (assuming that it is radiative) as well 
as in the g-wave propagation zone. In both zones, we have 
|Ci?l ~ ICvl S> max(|Cr|, ICpD) which justifies simplification of 
Eqs. (9) and (10). The relative scaling of ("r and is differ- 
ent in the evanescent and g-wave zone. In the former, where 
Vg and A are large, |Cr| ~ \Cp\ ~ sCi? ~ ^C<l' and the simpli- 
fication in Eq. (8) follows from j4 » 1. In the propagation 
zone, where the mode amplitude varies rapidly, wc have 



dr 



ICpI 



-ICr-l » ICvl > ICr 



which justifies the neglect of all terms containing S7. Note 
that the u)'^ term in the coefficient at could have been 
ignored too but traditionally this is not done. 

With the traditional approximation, the r and de- 
pendence in the eigenfunctions may be separated. The lat- 
ter dependence is given by the Hough functions. Following 
Townsend (2003a), wc substitute in Eqs. (8-11) 



C.=2yi(r)e(^), Cp = y2{r)ew 



(13) 



smi3Ci> ^ z{r)Q{-d), and isini?C^ = «(r)e(j?) (14) 
to obtain ordinary differential equations for yi and 3/2, 

.^ = (^,-3),,+ (^-V,),„ (15) 



and 

where, for a uniform rotation, A is the separation parameter. 
The radial dependence for the horizontal components of ^ 
is given by 

^_ 9 



(17) 



Eqs. (15) and (16) look just like in the case of no rotation, 
except that now we have w = Woha — mQ and in the place of 
£{£+ 1) the separation parameter A, which is determined as 
an eigenvalue in the equations for the Hough functions. We 
write these equations in the form 

e = -me + snQ (18) 
(l-/i^)^ = -msMe-|-(sV-l)© (19) 



{l-H^)^ = [X{l-iJ,'')-rn]e + msne, (20) 

where /i = cost? and s = 2^1 /uj is called the spin parameter. 
Eqs. (19-20) differ from Townsend's Eqs. (21-22) only in the 
sign of s, which he denoted v, in consequence of different 
sign convention for the time dependence. The equations to- 
gether with boundary conditions at /i = (symmetry) and 
/X = 1 (regularity) define the eigenvalue problem on A. The 
dependence of A on s and related properties of the Hough 
functions have been discussed in great detail by, for example, 
Bildstein et al.(1996), Lee & Saio (1997), and by Townsend 
(2003a). 

Of our interest here are only modes which are propagat- 
ing in radiative zones, where A > 0, that is, corresponding 



to A > 0. These arc the g modes for which A i{i + 1) at 
s — » 0, as well as the mixed gravity- Rossby modes, for which 
A changes sign from minus to plus at s = |m| -J- 1 (for brevity, 
in this paper, we shall call them r modes). It was shown, in- 
dependently by Townsend (2005b) and Savonijc (2005) and 
confirmed by Lee (2006), who did not use the traditional 
approximation, that such modes may be driven in B-type 
and early A-type stars. Let us note that if O is a function of 
r, then such modes may be trapped in the region of rapid 
rotation where A > while it is < in the rest of the star 
interior. 

In the case of massive main-sequence stars, the tra- 
ditional approximation cannot be used in their convoctivc 
cores. However, low frequency waves do not propagate there 
and thus the core surface may be treated as a boundary. Ex- 
cept for zero-age main-sequence stars there is a nearly dis- 
continuous transition at the core edge, r = rc, from N ^ 
to N ~ 0. In the core, i^p and other eigenfunctions arc slowly 
varying with distance from the center. Thus, with N ^ 0, 
Eq. (8) implies 



ICpI ~ e|Cr-| for r ^ rc 



(21) 



On the other hand, in the propagation zone, we have fc^Cp ~ 
N^Cr/g, where the local radial wave number is given by kr = 
kHN/uj and kn = v'A/r is the horizontal wave number. 
Hence, 

(22) 

Since there is no jump of density at r = rc, the continuity of 
5p and imply the continuity of ("p. Thus, from Eqs. (21) 
and (22), we get the boundary condition consistent with 
the traditional approximation to be imposed on solutions of 
Eqs. (15) and (16), 



2/2 (rc) 

?/l(fc) 



Ve» 0. 



(23) 



In stars of our interest, there are thin convective layers, con- 
nected with the Hell ionization and, if M > 7Mq , with the 
Fe opacity bump. Because of their small vertical extent, such 
layers cannot affect significantly global mode geometry. 

The two boundary conditions, which should be applied 
on solution of Eqs. (15-16), are 



2/2 = at r = rc and j/2 = j/i at r = R. 



(24) 



The same boundary conditions remain valid when non- 
adiabatic effects are taken into account. 

Eqs. (15-20) remain valid in the case of shellular rotation 
but then, both uo and A must be regarded as functions of r. 
The derivation given by Dziembowski & Kosovichev (1987) 
employed expansion of the 9 dependence in a series of the as- 
sociated Legendre functions. Assuming Q,/N ~ 5 -C 1, they 
decomposed the infinite system of differential equations in 
r into asymptotically independent second-order equations 
equivalent, in the limit 5 — + 0, to Eqs. (15-16). In their 
derivation there is an inconsequential omission of terms 
~ (5rflnft/dlnr and an implicit assumption that the log- 
arithmic derivative of ft is of the order of 1. 

Townsend (2005a) and Savonijc (2005) extended the 
traditional approximation to the case of non-adiabatic oscil- 
lations. Using the same equation as in the case of no rotation 
with only A in place of £{i + 1) is incorrect in the expression 
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for the horizontal heat losses. However, these losses, which 
are proportional to kjj, make rather small contribution to 
the overall work integral. Another simplification is the use of 
A with uj calculated in the adiabatic approximation. This is 
well justified in our applications because the non-adiabatic 
change in u) is very small. 

Lee & Saio (1989) pointed out that even at 5 <C 1 cou- 
pling between modes of close frequencies, the same m and 
symmetry may have a significant effect. In a specific ex- 
ample of 4AIq star sequence, Lee (2001) found that some 
modes, which are unstable if the traditional approximation 
is used, are found stable if the approximation is abandoned 
and the truncated expansion in Logendre functions is used 
for the eigenfunctions. However, many modes remained un- 
stable and we expect that this applies also in the case of 
models considered in this paper. We will discuss this matter 
further in Section 3. 



3 UNSTABLE MODES 

With the traditional approximation, stability properties of 
slow modes arc determined by two parameters: A and cu. 
In this section, we assume uniform rotation, hence lo is the 
mode eigenfrequency in the corotating system. The whole 
infiuence of rotation and the azimuthal order is absorbed 
in the A(s) dependence. Thus, like in the case of no rota- 
tion, we deal with a two-dimensional instability range for 
a specified model. In the present case, the two parameters 
are oj and s/X/oj. The latter determines radial order of the 
mode, n, and the r dependence of the Lagrangian pressure 
perturbation, Sp/p. For instability, the amplitude of the per- 
turbation should be large and slowly varying in the layer of 
the iron opacity bump. The frequency, uj, determines result- 
ing perturbation of the radiative flux. Furthermore, a match 
is needed between the thermal time-scale of this layer, rth, 
and the pulsation period, H = 2n/u). 

The instability ranges in terms of i and uj for non- 
rotating main-sequence B-stars were discussed in detail by 
Dziembowski, Moskalik & Pamyatnykh (1993) and by Pamy- 
atnykh (1999). A corresponding discussion for rotating stars 
was presented recently by Townsend (2005a). For a family 
of stellar models, VX/ijJ does not determine mode order but 
still, to a large extent, determines the shape of the pressure 
eigenfunction in the outer layers, which bring the main con- 
tribution to the work integral. Thus, the evolutionary radius 
increase results in longer periods of unstable modes. The as- 
sociated decrease of the effective temperature, causing an 
increase of rth, acts in the same direction. This is why the 
initial effect of evolution is an enhancement of instability. 
The tendency is reversed only near the end of the main- 
sequence evolution when dissipation in the g-wave propagar 
tion zone becomes significant. Too small values of rth are 
the cause why the modes with low values of A are stable in 
massive (M > A.SMq) hot objects. 

In this paper we present numerical results for a sin- 
gle stellar model calculated with Warsaw-New Jersey stellar 
evolutionary code (see e.g. Pamyatnykh 1999) assuming no 
mixing beyond the convective core and including only mean 
effect of the centrifugal force. The same approximations were 
adopted by Walker et al. (2005). Our model is characterized 
by the following parameters: M = 6Mq ,Xo =0.7, Z = 0.02, 



Xc = 0.307, logTcff = 4.226, Urot =300 km/s, which are very 
similar to those used by these authors. 

Our non-adiabatic pulsation code uses the same ap- 
proximation as adopted by Townsend (2005a) except that 
we apply the inner boundary condition 2/2 = at the edge 
of convective core. The code is almost a trivial modifica- 
tion of our standard nonradial pulsation code in its version 
adopting the Cowling approximation, which is indeed very 
accurate for all considered modes. In the adiabatic part of 
the code values of A arc interpolated from the tables of A(s) 
prepared with a separate code for all types of modes of our 
interest. The values were not changed in the non-adiabatic 
part. We focused the main attention on modes which may 
be driven by the opacity mechanism and are potentially de- 
tectable by photometry. Thus, the modes considered have 
A greater than zero but not too large. For each |m| ^ 2, 
we included g modes with the lowest values A. Sequences 
of modes of different radial order are denoted with {£,m), 
where (.{(.+1) = A(0). We also included two sequences of low- 
est order retrograde r modes with A > for s > \m\ + 1. The 
mode instability is characterized by the normalized work in- 
tegral 



where W is the usual work integral. The growth rate, ?5(ij.'), 
has the same sign but its value is affected mostly by mode 
inertia. The value of 77, which varies between -1 and 1, is a 
better measure of robustness of the instability and a better 
predictor of the mode amplitude. Fig. 1 shows the values of 
77 for all considered modes in their range of instability and 
close to it. Table 1 provides more information about most 
unstable modes (largest jy). 

The instability range in the observed frequencies, Vohe = 
Wobs/27r, extends from nearly zero for retrograde £ = 1 g 
modes to 4 c/d for prograde £ = 2 modes. The former modes 
have frequencies in the vicinity of Vrot = 0/27r. The modes 
with frequencies i/gtar = u!/2tt < z/rot are seen as prograde 
modes. The range of frequencies in the corotating system 
is much narrower. We may see in Table 1 that the most 
unstable modes occur near fstar ~ 1-6 c/d, A ~ 30 and 
n = 20. The prograde i = 1 modes have A < 2 and, thus, at 
VX/ijJ values favourable for driving, too long periods. Such 
modes were found unstable in somewhat cooler models. A 
lower value of Tgg would also result in lower i^star of the 
unstable modes. The shift in Vohs is the same but, if rnvrot > 
I'stB.T, of opposite sign. In the vicinity of our model a shift 
of 0.01 in Tgff leads to the 0.1 c/d frequency shift of most 
unstable modes. A change in the adopted value of Viot would 
shift the instability ranges in i^obs for m 7^ modes. 

There are many more strongly unstable modes than 
shown in Fig. 1. The instability continues up to A w 600. 
Modes with rj > 0.2 are found up to A « 120, which cor- 
responds to ?n = 11 for sectorial prograde modes. Owing 
to higher A values such modes are subject to much larger 
observable amplitude reduction. Visibility of modes is the 
subject of the next section. Here, we just want to point out 
that we should expect in any of multiperiodic B-type star 
many more excited modes than those detected. Such modes 
have properties similar to those listed in Table 1, except they 
would have much lower surface averaged amplitudes. Possi- 
ble excitation of many unobservable modes must be kept in 
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Figure 1. The instability parameter, rj (see Eq.(25)) as function of frequency for slow modes of low azimuthal orders in the model of 
HD 163868. The frequencies are calculated in the inertial system. The frequency of rotation is i/rot = 0/27r = 1.42 c/d. 



Table 1. Parameters of the most unstable modes with |m| ^ 2 for the stellar model with M = Q.OMq, log T^g = 4.226, log L/Lq 
and a rotational velocity of 300 km/s, which corresponds to Urot = 1.4165 c/d. See Eq.(29) for the definition of /. 
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mind while considering the feedback of oscillations on stellar 
rotation. The results shown in Fig. 1 and in Table 1 might 
suggest that there is preference to retrograde mode excita- 
tion but it is not found in the full set of the unstable modes. 
For instance, rj ^ 0.2 was found in over 450 progradc modes, 
and in less than 370 retrograde modes. Certainly, there is 
no evident asymmetry in driving prograde and retrograde 
modes. 

Our findings regarding mode instability difi^er from that 
presented by Walker et al. (2005). Unlike these authors, we 
do not find instability of the prograde i = 1 modes. On the 
other hand, we find many retrograde i = 1 modes, which 
may account for the longest periods measured in HD 163868, 
whereas no instability of such modes was found in that work. 
The general pattern of the difference is that we find insta^ 
bility in many more modes with higher A. Perhaps the tra- 



ditional approximation, which we rely on, is not applicable. 
The mode coupling, as described by Lee (2001), may sta- 
bilize only certain modes. Could it stabilize all retrograde 
i = 1 modes and destabilize all prograde i = 1 modes? We 
do not know the answers. Perhaps, there is a difference in 
internal structure between the models used by us and that 
used by Walker ct al. (2005). We should also be aware that 
the Lee & Saio method, which is based on a truncated series 
expansion, used in that paper, may not be accurate. 



Considering possible identification of modes detected in 
HD 163868, we should take into account visibility conditions 
for modes described by the Hough functions. 
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Figure 2. Light amplitudes in the Geneva V band (upper panels) and in radial velocity (lower panel) as functions of the inclination 
angle for modes listed in Table 1. 



4 VISIBILITY 

Observed amplitude of a mode depends on its true ampli- 
tude, geometry and the aspect angle, i. Photometric am- 
plitudes for modes described by single spherical harmonics 
were calculated in a number of papers. Results of corre- 
sponding calculations for modes described by the Hough 
functions were first published by Townsend (2003b), who 
relied on truncated expansion of Q in associated Legendre 
functions. In our calculations of the photometric and radial 
velocity amplitudes, we avoid such expansion and carry out 
two-dimensional integration over visible hemisphere. Details 
of our approach and a survey of the results will be published 
elsewhere (Daszyriska-Daszkiewicz et al., in preparation). 
Here we quote only the background formulae and present 
discussion of visibility of modes listed in Table 1 as a func- 
tion of aspect. 

Within the linear non-adiabatic theory, there is a free- 
dom in normalization of eigenfunctions. Here we adopt the 
most common one, that is yi{R) = 1, and 



e^d/i = 1. 



Thus, we write 



^r{R) = eROZ, 



(26) 



where Z = exp[i{mip — u}t)]. The value of |e| (not to be con- 
fused with e in Section 2) determines true mode amplitude. 
Note that the surface displacement is given by 5R[^r(i?)] and 
its rms value is equal to |£|/-\/47r. The surface boundary con- 
dition, y2{R) = yi{R), combined with Eqs. (4) and (14) leads 
to the following expressions for the horizontal components 



of the displacement 



and 



_ GM e 

' uj'^R'^ sint 



GM e 



(27) 



(28) 



uj^R'^ sin 1? 

These components are needed only for calculation of the 
pulsation velocity field. 

For evaluation of the photometric amplitudes, we need 
a linear relation between the radial components of displace- 
ment and perturbed flux which we write in the form 



efQZ, 



(29) 



where / is a complex coefficient determined from linear non- 
adiabatic calculations. Its values for modes considered arc 
given in the Icist column of Table 1. In order to calculate the 
flux in a specifled band, x, we need also an expression for 
the perturbed surface gravity, which is 

5g _ f ruP 
9 \ 9 

and the perturbed surface element. 



2 ez w 2ez, 



5dS 

dS 



29, 



ae ime 



where dS = R^djidifi. The tangential components arise from 
the change of the normal to the stellar surface. These com- 
ponents are needed also in the expression for the perturbed 
limb-darkening, 5hx- The perturbed flux is thus given by 



5L, 



Js 



hx + Tx5hx)dSer + Pxhx^dS] ■ nobs, (30) 
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where X represents a photometric band, riobs is the unit 
vector toward observer, and the integration is done over the 
visible unperturbed hemisphere. The hmb-darkening coeffi- 
cient is normalized, so that 

Js 

where /i = nobs • e,-. Static atmosphere models are used 
to calculate wavelength-dependent flux, Tx{Tet{,g), and the 
limb-darkening coefficient, hx(Tcit,g). 

The disc-averaged change of radial velocity is given by 
the expression 

<'yrad >= / (Wpuls + Vrot) ' Hobs^x/lx Hobs • dS, (31) 

Js 

where 



= — and 
dt 



Vrot 



= RQ sin i?eu: 



In the linear approximation, the contribution from the first 
term may be calculated with unperturbed h, J^, and S. The 
contribution from the second term arises solely from pertur- 
bation of these quantities. Thus, the radial velocity change 
is given by 



<'yrad> 



-a 



&Tx , 



Vvot 



(Wrot • "^obs)-^^ • Hobs }J^xhxdS. 



■ Hobs A + 

(32) 



We show in Fig. 2 the amplitudes of light and radial 

velocity as functions of the inclination angle for the modes 
listed in Table 1. These amplitudes are calculated with the 
arbitrary choice |e| = 0.01 and for the Geneva V band, which 
is similar to Johnson's V . The plotted quantities are 



Av = 1.086 



5Lx 



and 



-^vrad — I *^ ^rad ^ 



In the dependence Avii), we see certain patterns which 
are the same as for modes described by single spherical har- 
monics. Zonal modes are best seen from near polar direction, 
while symmetric (even) modes are best seen from near equa- 
torial direction. There are, however, significant differences. 
In particular, an increase of A does not lead to so large am- 
plitude reduction as expected from the — > A replace- 
ment. Compare, for instance, two sectorial i = 1 modes. 
There is the factor 25 difference in A between m = — 1 
and m = +1 modes and only the factor 4 in Ay, while 
the factor ~ 12 would be expected in the non-rotating 
case, on the basis of the decline of the disc-averaging fac- 
tor, bt, with I (e.g. Dziembowski 1977). According to our 
calculations the m = -1-1 mode is stable. Among unstable 
modes, the {i = l,m = —1) is the second best visible, af- 
ter {£ = 2,m = +2), from the near equatorial directions. 
Despite of near-equatorial trapping caused by rotation, the 
optimum inclination for detection of tcssoral modes, such 
as (^ = 2,m = ±1), is shifted toward the pole because 
then cancellation of contributions from the two hemispheres, 
which have opposite signs, is reduced. Such a cancellation is 
responsible for rather low amplitudes of the r modes. 

The behaviour of the radial velocity amplitudes, de- 
picted in the lower panels of Fig. 2, shows rather different 
pattern. At intermediate inclinations of the rotation axis. 
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Figure 3. Growth rates for symmetric (even) modes, which are 
preferentially visible from the near equatorial directions (see Fig. 
2) , in models of HD 163868 calculated with two indicated equato- 
rial velocities of rotation (the upper panels) and amplitudes the 
peaJis detected in the MOST data. 



the r modes turn out most easily detectable. At large inclina- 
tions (near equatorial observer) the modes \m\ = £ are best 
seen both in light and radial velocity variations. However, 
the ratio exhibits a strong mode dependence. Interpreta- 
tion of the radial velocity amplitude pattern is complicated 
because the contribution from pulsation and rotation are 
of comparable size. Here we would like only to stress that 
the strong mode dependence of the amplitude ratios points 
to the good prospect for mode identification by combining 
spectroscopy and photometry data. 



5 PEAKS IN THE HD 163868 OSCILLATION 
SPECTRUM 

The measured value of v sin i in HD 163868 is high. There- 
fore, we assume that the star is seen from a near equatorial 
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direction and, as possible candidates for identification, we 
consider only symmetric modes. 

There are three groups of frequencies found by Walker 
et al. (2005) in the star. That with Fobs < 0.5 c/d could 
only be interpreted in terms of the retrograde £ = 1 modes. 
Possible interpretations for the group with Vohs between 1.4 
and 2 c/d are £ = 2 modes with m = —2 and m = 0. For the 
group with fobs between 2.8 and 3.75 the only interpretation 
is £ = 2,m = +2. In this last case our mode identification is 
the same as that of Walker ct al. (2005). 

In Fig. 3 we see that our interpretation is also not fully 
satisfactory though, unlike our predecessors, we have no dif- 
ficulty with interpretation of lowest frequency group, where 
the highest amplitudes arc found. .Just like theirs, our pre- 
dicted range at the highest frequencies is somewhat narrow. 
The largest discrepancy between results of our calculation 
and observation is for the intermediate group. Certainly, the 
discrepancy may be reduced by adjusting Vrot and Toft. 

It is possible that the some of the peaks in the HD 
163868 osillation spectrum, in particular those in the low 
frequency range, are artefacts of the data processing. In- 
terpretation of peaks in terms of specific modes should be 
tested with simultaneous multiband photometry and radial 
velocity data. We have seen in the previous section that the 
ratio of radial velocity to light amplitude is strongly mode 
dependent and it is exceptionally large for the retrograde 
i = 1 modes, which may explain peaks in the low frequency 
range. The dominant peak in this range and in the whole os- 
cillation spectrum has photometric amplitude of 0.013 mag 
in the MOST photometric system. Assuming that the am- 
plitude in the V band is similar and using plots in Fig. 2, 
we find radial velocity amplitude of about 30 km/s, which 
should be easily mcasm-ed in radial velocity data spanning 
sufficiently long time interval. 



6 CONCLUSIONS 

Data on oscillations in Be stars may be helpful in explaining 
activity of these stars. The application is twofold. Firstly, the 
data may be a source of seismic information on internal ro- 
tation. Secondly, the data may allow to Eissess the role of the 
feedback effect of oscillation in these stars. Both aims require 
identification of modes detected in the observed oscillation 
spectra. The clue to identification which may be provided 
by the theory is the determination of a set of modes that 
may be excited and detected. 

Effects of rotation, which are important in Be star oscil- 
lations, must be taken into account. This is most easily done 
with the traditional approximation. Validity of this approx- 
imation in application to slow modes in Be stars, we dis- 
cussed in Section 2. We applied our non-adiabatic version 
of this approximation to a model of HD 163868, the first 
Be star with the rich oscillation spectrum. The spectrum, 
determined from the MOST photometric satellite data, and 
its preliminary partial interpretation, has been recently pub- 
lished by Walker et al. (2005). 

We showed that there is a large number of unstable 
modes in the model of this star, covering the whole frequency 
range of the peaks in the HD 163868 spectrum and beyond. 
The instability extends to modes of high horizontal degrees, 
corresponding to £ w 24 in the case of no rotation. We found 



no significant difference in driving between prograde and 
retrograde modes. 

Our calculations of the relative visibility of the unstable 
modes show that the three consecutive groups of the peaks 
may be explained by low-degree high-order g modes seen 
from the near equatorial direction. Specifically, the group 
with frequencies below 0.5 c/d may be explained by retro- 
grade £ = 1 modes, the group axound 1.7 c/d by retrograde 
sectorial and zonal £ = 2 modes the group around 3.5 c/d 
by prograde sectorial £ = 2 modes. Only the interpretation 
of the highest frequency group is the same as that proposed 
by Walker et al. They interpret the intermediate group in 
terms of £ = 1 prograde mode, which we find stable. In 
contrast, the modes that we associate with this group are 
found stable in their calculations. Also they found all retro- 
grade £ = 1 modes to be stable and therefore did not find 
a satisfactory interpretation for the lowest frequency group, 
where peaks have the highest amplitudes. They proposed 
that some of the peaks in this group may be interpreted in 
terms of retrograde r modes, which are found unstable, both 
by them and us, but their interpretation leaves unexplained 
highest peaks and there is also a problem with the visibility 
of such modes from the nearly equatorial direction. A possi- 
ble cause of the difference in interpretation of the same data 
is the use of very different approximations in the treatment 
of the difficult problem of slow oscillations in rotating stars. 
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